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Abstract. The global properties of static perfect-fluid cylinders and their 
external Levi-Civita fields are studied both analytically and numerically. The 
existence and uniqueness of global solutions is demonstrated for a fairly general 
equation of state of the fluid. In the case of a fluid admitting a non-vanishing 
density for zero pressure, it is shown that the cylinder's radius has to be finite. 
For incompressible fluid, the field equations are solved analytically for nearly 
Newtonian cylinders and numerically in fully relativistic situations. Various 
physical quantities such as proper and circumferential radii, external conicity 
parameter and masses per unit proper/coordinate length are exhibited graphically. 
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1. Introduction and summary 

Since 1917 when Levi-Civita gave his static vacuum solution, cylindrically symmetric 
spacetimes have played an important role in general relativity. Recently attention has 
been mainly paid to dynamical situations or quantum issues. Some of these aspects 
are surveyed in [1] (section 9 and references [247-263], [275-277] therein), for example; 
for even more recent works, see, e.g., [2, 3, 4]. 

Still, in a recent concise review on static cylinders [5], Bonnor introduces the 
topic by stating that 'the Levi-Civita spacetime continues to puzzle relativists'. In 
contrast to the Schwarzschild metric described completely by a single parameter, the 
Levi-Civita solution contains two essential constant parameters — m, related to the 
local curvature, and C, determining the conicity of the spacetime. Relating these two 
parameters to physical sources turns out to be delicate, as we recently found in the case 
of static cylindrical shells of various types of matter [6] , and as we will demonstrate 
for solid perfect-fluid cylinders in the present paper. 

In the second edition of the 'exact-solutions bible' [7], there are quoted about 20 
papers on static perfect-fluid cylindrically symmetric fields. However, these solutions 
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are local, no analysis of global properties is usually available. In addition, most perfect- 
fluid equations of state in the literature are either ad hoc or not very physically 
plausible. In fact, the cylindrical counterpart of the famous Schwarzschild interior 
solution, in which the fluid is incompressible, so that its matter density \x = = 
constant, is not known. Although this equation of state is unphysical (implying infinite 
velocity of sound), in the spherical case not only can the corresponding solution be 
easily found but it also gives 'not unreasonable' estimates for maximum masses of 
neutron stars, for example. 

The purpose of the present work is to study static perfect-fluid cylinders and 
their fields globally, to investigate the Newtonian limit and to analyse incompressible 
cylinders numerically. In section 2, we start from a line element in which the 3-space 
is rescaled by the norm of the timelike Killing vector and we choose the distance 
from the axis as the radial coordinate. Usually we require the axis to be regular, 
but the presence of a conical singularity (a 'cosmic string') along the axis is also 
considered. The field equations in these coordinates enable one to find convenient 
relations between some of the metric functions. In particular, the vacuum Levi-Civita 
solution can be rederived in these coordinates. We show in detail that although it 
contains five constants of integration, only two — the mass parameter and the conicity 
parameter — cannot be removed. 

In section 3 we first formulate a lemma (and then prove it in appendix A) stating 
that for a smooth equation of state ji{p) and any value of the central density /io there 
exists a unique solution of the field equations in a neighbourhood of the axis, which 
can be continued into a global solution. Depending on the equation of state, the fluid 
occupies the entire space, or the cylinder is of a finite extent in the radial direction 
provided that the pressure vanishes at a finite radius, p{R) — 0. Another lemma 
shows that a unique vacuum solution can be joined to the inner perfect-fluid solution 
at this radius R where the pressure vanishes. Finally, a Theorem is proven asserting 
that systems with fi(p = 0) > always have a finite radius — like with spherically 
symmetric balls. 

The field equations are written in terms of A = 1/c 2 (c — the velocity of light); 
for A = the Newtonian limit is recovered. The conformal 3-space becomes flat 
and the relativistic Tolman mass turns into the Newtonian mass per unit length. 
This is demonstrated in section 4. In appendix B, the Newtonian static constant- 
density cylinders and spheres are constructed directly within the Newtonian theory 
for comparison. (For shell sources in Newtonian gravity, see appendix in [6].) 

In the brief section 5 it is explicitly shown how an outer Levi-Civita solution 
can be joined to an inner perfect-fluid solution at the boundary where the pressure 
vanishes. The need of a non-trivial conicity parameter, C / 1, outside a general 
perfect-fluid cylinder is discussed in detail in appendix C. As an example, the Evans 
solution [8] with the equation of state fi = fio + 5p, > 0, is employed to elucidate 
that only special cylinders can be joined to the outer Levi-Civita metrics with the 
conicity parameter C — 1. 

Section 6 is devoted to cylinders of incompressible fluid. As mentioned above, 
exact solutions representing incompressible perfect fluid cylinders have not yet been 
found. We first study analytically weakly gravitating cylinders. Expanding all 
functions entering the field equations in terms of the dimensionless radial distance 
from the axis, the pressure at any point is found in terms of the central pressure 
p c . Inverting the series, the radius of a cylinder can be determined by requiring the 
pressure to vanish. In this way, relativistic parameters characterizing the cylinder 
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can be determined analytically as corrections to the Newtonian values. The conicity 
parameter starts to deviate from its Minkowskian value only by terms of the order 
0(Ul lnll c ) where II C = p c /fi c 2 . 

Fully relativistic cylinders of incompressible fluid are treated numerically. Here 
we extend considerably the results of Stela and Kramer [9]. We allow the conicity 
parameter C ^ 1 (as one should — see above) , admit high values of the central pressure 
and we also exhibit other quantities of interest graphically. In particular, we show how 
the proper radius of the cylinder increases with increasing central pressure, whereas the 
circumferential radius starts to decrease. The peculiar behaviour of the circumferential 
radius is well seen in the embedding diagram of the (z, t) = constant surfaces and their 
dependence on the central pressure. Both the external Levi-Civita mass parameter 
and external conicity parameter are plotted as functions of the central pressure. It 
is illustrated graphically that the dimensionless mass per unit proper length has an 
upper bound of 1/4, whereas Thome's C-energy scalar is bounded from above by 1/8. 
We also plot dimensionless mass per unit coordinate/proper length of the cylinders as 
functions of their proper and circumferential radii and compare them with the case of 
incompressible spherical stars. 

We believe that full understanding of static cylindrical systems will give intuition 
for more realistic situations in the neighbourhood of elongated but finite systems. In 
addition, the only realistic hope to construct exact examples of interaction between 
simple physical matter and gravitational waves seems to live in cylindrical symmetry. 
Even perturbative calculations of slowly rotating and oscillating cylinders might bring 
useful insights. The simple static configurations would then serve as convenient 
backgrounds. 

2. Field equations and their vacuum solutions 

We define cylindrical symmetry for static non- vacuum spacetimes by assuming that in 
addition to the timelikc Killing vector, two further commuting Killing vectors acting in 
hypersurfaces orthogonal to the 'static' Killing vector exist 1 . Adapting coordinates to 
the symmetries, we can write the line element with the metric in the 3-space rescaled 
by the norm of the timelike Killing vector (written as an exponential) in the form 

ds 2 = -e^c 2 dt 2 + e-^(A 2 dr 2 +B 2 dv 2 + C 2 dz 2 ) , (2.1) 

with A(r), U(r), B(r), C(r), where we keep the velocity of light c because we wish to 
discuss the Newtonian limit in detail. We assume circular orbits for Jj^ and consider 
only spacetimes with a regular axis. Occasionally, we admit an infinitely thin cosmic 
string, which is described by a conical singularity along the axis. The distance from 
the axis (in the conformal 3-metric) can be used to define a unique radial coordinate, 
i.e., we put A = 1. There are many other geometrical choices for the radial coordinate. 
If we fix the range of ip to be [0, 2tt), we obtain the following regularity conditions on 
the axis by comparing our form of the metric with the Minkowskian metric (possibly 
containing the conical singularity) written in cylindrical coordinates: 

1/(0) = 0, U'(0) = 0, 

5(0) = 0, B'(0) = l/C*, (2.2) 

C(0) = 1, C'(0) = 0, 

1 For the definition of cylindrical symmetry, see [7]; for a detailed, careful discussion, see, e.g., the 
recent work [10]. 
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where C* is the axis-conicity; in particular, the axis is regular for C* = 1. The functions 
U,r~ 1 B and C have to be diffcrcntiable functions of r 2 to ensure differentiability on 
the axis. The values (2.2) of U and C on the axis determine the normalization of the 
static and translational Killing vectors. 
For fluid cylinders we have 

T afj = (c 2 fi + p)u a Uf3 + V9 a fi , (2.3) 

where the 4- velocity u a is normalized as 

g a pu a u fi = -1. (2.4) 

The fluid is assumed to be at rest in the coordinates of metric (2.1), hence u a = 
(c _1 exp(— U/c 2 ), 0, 0, 0). The Einstein field equations 

8vrG. 1 . 8nG . . 

Roc/3 = —^{T a p - -Tg a p) <=> G a p = —^T a p , (2.5) 

imply (A - 4,) 2 : 

Ru = U"BC + B'CU' + C'BU' = AitG{/i + X3p) e - 2UX BC , (2.6) 

G rr = -X 2 (U') 2 BC + B'C' = 8nGe- 2UX BCX 2 p , (2.7) 

G w = \ 2 (U') 2 C + C" = 8irGe- 2UX C\ 2 p , (2.8) 

G zz = \ 2 {U') 2 B + B" = 8irGe- 2UX B\ 2 p . (2.9) 

As we have three unknown functions and three second order equations (2.6), (2.8), 
(2.9) we expect six integration constants. The first-order equation (2.7) reduces the 
number of constants in the solution to five. The conservation law T a f$ ;a =0 gives 

U' = -t^— ■ (2.10) 
Xp + ii y ' 

We assume that an equation of state (EOS) is given. Then we can integrate the 
last equation, 

UU-UW-UM-fj^. (2.1!) 

For incompressible fluid of density ^ an d central pressure p c , this can be integrated 
explicitly to yield 

U(R) = Iln(l + A^) (2.12) 
A 

for the value of U on the surface r — R with p(R) = 0. The field equations can be 
simplified as follows: from equations (2.7), (2.8) and (2.9) we can eliminate U' to 
obtain 

= 8nGX 2 e- 2U Ap. (2.13) 

BG 

Equation (2.6) can be rewritten as 

U" + ^^U' = 4irG(n + \3p)e- 2UX . (2.14) 
BC 

2 To avoid confusion, hereafter we denote the inverse square of the velocity of light by A. Since C is 
a metric function, we denote the 'external' Levi-Civita conicity parameter by C (cf. equation (2.19) 
and below) and the 'internal' axis-conicity by C* as in equation (2.2). 
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Hence, we obtain two equations for U and for the product BC. 

One immediate consequence of the field equations: take Cx (2.9)— Bx (2.8) and 
consider the above initial conditions (2.2). This yields 

CB BC ' = ^ = constant (2.15) 

for < r < co. 

To construct local solutions having a prescribed EOS and ignoring the regularity 
of the axis, one can proceed as follows: given p,(p), equation (2.11) implies there are 
unique functions /x({7), p(U). Inserting these functions in equations (2.13), (2.14), 
we can solve for U and BC locally. Next, equations (2.8) and (2.9) determine B 
and C. Equation (2.7) is satisfied as a consequence of the Bianchi identities. Hence, 
given an EOS, there are many local solutions away from the axis. Infinitely extended 
solutions with some special equations of state (fj, = -fp) were derived in [11, 12, 13], 
while finite solutions satisfying energy conditions were obtained in [14] and [8]. In 
[15], the authors consider polytropic equations of state using numerical solutions of 
the structure equations. 

We can also easily find vacuum solutions in the above coordinates. For simplicity, 
we now put A = 1/c 2 = 1. Equation (2.13) implies BC = ai(r + a 2 ) with constant 
parameters a\ ^ and a 2 (if a\ = then either B = or C = in contradiction with 
(2.1)). From (2.6), we have {U 1 BC)' = 0, so we can write U'BC = 0103 where we 
introduced a third integration constant CI3 3 . Substituting for BC we find 

U-«*(!±±). (2.16) 

with 04 7^ being a fourth integration constant. To find B and C, we have to solve 
(2.8) and (2.9) with vanishing rhs, taking into account BC = ai(r + a 2 ). The solutions 
can be written as 

where 05 7^ is a fifth integration constant and n^ 2 = (1 ± \J\ — 4o|)/2. (We can 
also interchange the roles of B and C . As long as we have no axis we cannot really 
distinguish between the two Killing vectors d z and d v .) From the relation for ni i2 we 
can see that the range of a 3 for which solutions exist is 03 € [—1/2, 1/2], It is thus 
advantageous to introduce a new constant 7 € [0, 1] such that 03 — — 7). With 

this definition, the exponents in (2.17) can be written n\ = 7, n 2 = 1 — 7. Introducing 
5 = a\a,5, a = 0,5/0,4, b = a 2 and L = a 5 , the resulting metric functions read 

r + b\ „ / r + b\ 7 „ ^fr + b 



[/ = V / 7(r 3 7)ln U-h . C= -H > S = M-H • (2-18) 



The solution thus has five constants of integration. However, we now show that 
three of these constants are redundant. We shall see that this is related to the fact 
that the Killing vectors dt,d z are defined only up to constant multiplicative factors 
and, in addition, the whole metric can be rescaled by a constant. First, denote 

a = [a(l — and introduce new constants p , m and C by the relations 

L = a Po , 7=T^, 5=^(l + m 2 )^. (2.19) 
1 + m z C 

3 Note its relation to the Tolman mass, mj< of (4.8), following from U' BC = 2mj = constant. 
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Then apply the transformation 

2 / \ 1 + TO 2 

[t,z,r, a (l + m 2 )Tf^C, J^{fj - b ^l (2- 2 °) 

This leads to the metric of the form 

(\ 2m / \ 2m(m— 1) / \ —2m 

~ ) dr2 + v ) ^ + d ^ + 75 M ~ ) ^ ^ 2 - 21 ) 

containing three constants only. 

Next, use a linear coordinate transformation with a scaling parameter q, which 
will be fixed subsequently, as follows: 

m m(m-l) m(m-l) 

[t,C,P,v] - [g^^t, g^^z, g^^TTr, p]. (2.22) 

This yields 

/rn\ 2m /toX 2 ™^- 1 ) /i __™2_\ 2 / rn \ - 2m 

d -' - - 6) di2 + 6) " L - 2 + * 2) + (? G) <2 - 23) 

By choosing q appropriately, we can dispose of cither C or p so that the resulting 
metric involves only two constants. Indeed, setting q = p we eliminate p in (2.23): 

ds 2 = -r 2m dt 2 + r 2m ^ m - x \dz 2 + dr 2 ) + ^^dp 2 , (2.24) 

m 2 

where we put C — C / p^^ m+1 . With C = 1, this is the standard Levi-Civita metric as 
given in the primary reference [7], equation (22.7). In (2.24), the t, z and r coordinates 
and the conicity parameter C do not have the usual dimensions. We can maintain the 

conventional dimensions if we set q — C m 1 , thus obtaining 

/ r \ 2m / t \ 2m(m- 1) / t \ —2m 

ds 2 = -(£) dt 2 + (£j {dz 2 +dr 2 )+ r 2 ^ dp 2 , (2.25) 

~ m 2 -m+l 

with R = po/C m 1 . Notice, however, that this cannot be done for m = 0. Hence, 
this coordinate system has the usual dimensions but (2.25) does not include the 
standard cosmic string as the limit m — ► only yields the flat spacetime without 
deficit angle. In the following, we use the form (2.24). 

Note that if the symmetry is considered locally only, we can also rescale <p and 
thus put C = 1 in (2.24). This is not the case, however, as we insist on a particular 
range of tp and require <p £ [0, 2ir). 

3. Existence of global solutions 

Our aim here is to show that for a given barotropic EOS, a unique global solution 
corresponds to each value of the pressure on the axis p c > 4 . Depending on the EOS, 
the fluid cylinder either is of a finite extent, or the fluid occupies the entire space. (In 
this section we again put A = 1/c 2 = 1.) 

We begin with the construction of solutions with a regular axis. Inserting Taylor 
series for the metric coefficients into the field equations indicates that for a given 

4 We choose the pressure rather than the energy density as the fundamental variable since in section 
6 we analyze cylinders of incompressible fluids. 
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analytic equation of state n(p), a value p c of the pressure at the center determines a 
unique solution. We prove the following lemma in Appendix A. 

Lemma 1: Let fi(p) be a barotropic equation of state. Then for each value of the 
central density /i = m(Pc) > there is r > such that a unique solution of the field 
equations (2.6)-(2.10) exists for < r < r . 

Next we show that the solution can be extended as long as the pressure is positive. 

Lemma 2: Suppose we have a solution with a regular axis. Then U , B and C are 
monotonically increasing and p is monotonically decreasing as long as p > (then 
also [i + p > 0). 

Proof: B' and C are positive near the center. If any of them vanished at some 
point, then equation (2.7) would give a contradiction because we assumed p > 0, and 
B, C > near the center as a consequence of the regularity conditions (2.2). Hence, 
B and C are monotonically increasing. 

Equation (2.14) implies— due to U'(0) = 0— that U"(0) > 0. Hence U' is positive 
near the center. Assume U'(r ) — 0; then U"(r ) > means that U cannot have a 
maximum at ro, so U is monotonically increasing, i.e., U' > 0. Equation (2.10) then 
implies p' < 0. (This is really the maximum principle for AU — [i + 3p.) 

Now we have the following possibilities: 

(1) The pressure is positive for all values of r. Then the fluid fills all the space and 
it is not possible to extend the solution because radial (r-)lines have an infinite 
proper length (cf equations (2.1), (2.11)). 

(2) The pressure vanishes at some value of r. Then we show in the next lemma that 
a unique vacuum solution can be joined to the inner solid solution. 

Lemma 3: If p{R) ~ 0, a unique vacuum solution — a particular Levi-Civita 
solution — can be joined to the fluid cylinder. 

Proof: We first show that U(R) is finite. For equations of state with a positive 
boundary density, fj,b = fJ,(p = 0), this is obvious from equation (2.11). If however, the 
boundary density vanishes, the integral in equation (2.11) could diverge for p — > 0. 
If this happens, we use the fact that U,B,C, U', and (BC)' are positive. Therefore, 
from equation (2.14) we have 



This shows that U' is bounded for < r < R if the EOS is monotonic (in fact, 
it is sufficient to assume that /i(p) is bounded on [0,p c ]). The same holds for U. 
Furthermore, the monotonicity of U and the upper bound on U" imply that U, U' 
have limits at R. As U, U' have limits at R, we can consider equations (2.8) and (2.9) 
as linear equations for B and C with given coefficients. This implies that B, B' , C 
and C have limits at R. 

All field equations remain meaningful for (i = p = 0. Therefore, the constants 
B(R), B'(R), C(R), C'(R), U(R), U'(R) provide the boundary values for the vacuum 
field equations; these have a unique solution — a Levi-Civita solution. Since the metric 



U" < 4ttG(m + 3p) 



(3.1) 



so that 




(3.2) 
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is C 1 in our radial coordinate, the junction conditions at the boundary are satisfied 
by construction. 

Summarizing, we have shown that a barotropic EOS determines a 1-parameter 
family of global spacetimes. 

As mentioned in the introduction, solutions with constant density, the analogues 
of the interior Schwarzschild solutions, have so far not been found as explicit exact 
solutions. In section 6 examples of cylinders with ^=constant are constructed 
numerically. The following theorem shows that a general constant-density cylinder 
always has a finite radius. 

Theorem: For a solution with a regular axis, a barotropic EOS and a positive 
boundary density fi(p — 0) = > 0, there is a finite radius R with p{R) = 0. 

Proof: Suppose that p is positive for all r. Equation (2.11) implies that U has a 
limit for r — > oo. We want to show first that U' — ► 0. If this is not the case, then 
there is a sequence r n — ^ oo with U'(r n ) = a > 0. Since U' is integrable, the peaks 
where U' reaches the value a repeatedly must become ever narrower with increasing 
r's. Therefore, there must be values of r where there are arbitrary high positive and 
negative values of U" . However, equation (2.14) shows that U" is bounded and we 
have a contradiction. Thus, U' — > 0. 
We can rewrite equation (2.6) as 

(BCU'Y = ^G{p + 3p)e- 2U BC , (3.3) 

so that by integration we get 

[/' = J- / AirG(p + 3p)e- 2U BCdr . (3.4) 

Eliminating the pressure through equation (2.13), we obtain 

u ' {r) = icf ^ G ^ 2UBCdr + lkJo h BC Y' dr ' < 3 - 5 ) 

or 

U'(r) = ^f o ^G^BCdr+^l{BCy{r)-^l{BC)'{0) = 

= Jo ^G»e-™BCdr + ^|((BC)'(r) - 1/C*). 
Due to equation (2.15), we further have 

(BC)' - ^ =CB' + BC - ^ = CB' - BC + 2BC - ^ = 2BC' > 0. (3.7) 

If BC is bounded, the first integral gives a positive contribution to V and we have a 
contradiction with U' — > because the second term in (3.5) is also positive. For BC 
unbounded, we calculate the limit of the integral for r — > oo by l'Hospital rule as 

,. ^G[ie- W BC 

^ (Bey • (3 ' 8) 

Now for r — > oo either both terms in equation (3.6) (which are positive) are finite, or, 
if one vanishes, the other diverges and we again have a contradiction with U' — > 0. 
Hence, p has to vanish at some finite radius. Note that equations (3.6), (3.7) and 
(2.15) imply that this conclusion is valid for both a regular axis (C* — 1) and for an 
axis with a cosmic string (C* ^ 1). 



(3.6) 
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4. Newtonian limit 

For a systematic treatment of the Newtonian limit, we refer to [16, 17]. Above, we 
have already chosen our variables in such a way that the field equations have a limit 
for a diverging velocity of light, i.e., for A — > 0. For A = 0, equations (2.6-2.9) imply 

U"BC + B'CU' + C'BU' = A-kGiiBC , (4.1) 

B'C = , (4.2) 

C" = , (4.3) 

B" = , (4.4) 

U' = -- . (4.5) 

Equations (4.2)-(4.4) show that the conformal 3-space is flat in the Newtonian limit. 
To adapt the coordinates to the axial symmetry and the regularity of the axis (cf 
equations (2.2)), we have to choose B = r, C = 1. Equation (4.1) then becomes the 
Poisson equation for the Newtonian gravitational potential: 

U" + -U' = 4mG[i. (4.6) 
r 

The field equations depend on A in such a way that it is easy to demonstrate 
the existence of families of exact solutions that have a Newtonian limit. Suppose we 
choose an EOS that is independent of A. The theorem in Appendix A, which is used to 
obtain solutions with a regular axis, also shows that these solutions depend smoothly 
on A and that the limiting solution satisfies the equations for A = 0. The same is true 
for the extension of the solution up to the radius where the pressure vanishes. 

Suppose we have such a A-family of solutions with finite radii R\ for all A's. 
Integration of equation (2.6) up to the boundary gives 

(U'BC)(R X ) = G4tt [ V + \3p)BCe- 2XU dr . (4.7) 
Jo 

The constant 

m T = 2tt [ \p + X3p)BCe- 2XU dr (4.8) 
Jo 

is the Tolman mass (see, e.g., [5, 18]). Since all metric coefficients have limits for 
A — > 0, equation (4.8) implies that the limit of tut is the Newtonian mass per unit 
length 

lim ttit = 2ir nrdr . (4.9) 
*-»o Jo 

The corresponding results for Newtonian cylinders and spheres are summarized in 
Appendix B. 
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5. Joining the inner fluid solution and the outer Levi-Civita solution 

Next, we want to demonstrate explicitly how a fluid solution, with the property that 
the pressure vanishes at some fixed r — R such that all metric coefficients and their 
first derivatives are finite there, can be uniquely matched to a vacuum solution given 
by functions (2.18) 5 . 

This is obvious because all the equations are regular for p = 0. Hence we can 
just take the values of the metric and their first radial derivatives as the initial values 
for the vacuum field equations and obtain a unique solution satisfying the matching 
conditions (the metric and its first derivatives continuous) by construction. Let us note 
that, as we have checked, we obtain the same results by requiring that (i) the cylinder's 
proper circumference is the same as measured from both sides of the surface and (ii) 
the surface energy-momentum tensor (calculated using the general Israel formalism 
[19]) vanishes. 

Equations (2.15) and (2.18) imply 

On the surface, we have to satisfy equations (2.18): on the left-hand sides we have 
the metric potentials obtained by integration within the cylinder, from these we now 
determine the five constants on the right-hand side. 

Denoting F(r) = B(r)C(r), we have on the surface F'(R) = S/L = 1/C*(1 - 2 7 ) 
and thus 7 = (C*F ' - 1)/{2C*F ') = m 2 /{l + m 2 ). Using this, we readily find 



C* F ' — 1 

m = Vc^TT' (5 ' 2) 

where F' is evaluated at R. It is elementary to show that the outer Levi-Civita 
parameter m is bounded by to <G [0, 1). Using the expression for the Tolman mass 
m T = 2tt f(f e-^BC'ifi + 3p)dr, we find 

Thus, we obtain an explicit relation between the Levi-Civita parameter, to, of the 
outer metric, the axis conicity, C* , and the Tolman mass, tot, of the fluid cylinder: 

mT= 2^T^- (5 ' 4) 

Using the external metric in the form (2.24), we find the conicity parameter from the 
continuity of the metric and its radial derivative. Simple calculations yield 

C f' 



1 / B \ V(e*F') 2 -l-2C*F' 

c = v{b^} (5 ' 5) 



where B, B' , U and F' are evaluated at R, the axis conicity C* is constant. 
Consequently, the conicity parameter, C, of the outer metric is determined by the 
values of the metric potentials on the surface independently of the EOS. In the case 
of incompressible fluid, the coefficient e u on the surface is given explicitly by (2.12). 
Since the surface values of the metric, for a given EOS, are determined by the central 

5 In the following, we need to distinguish between L of (2.18), po of (2-21), R of (2.25) and the 
coordinate R, proper R p and circumferential R c radii of the cylinder defined by equations (6.8), 
(6.9). In this section again A = 1/c 2 = 1. 
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pressure p c , both m and C are, in fact, given by p c and by C* in case of a general 
cylinder. Even if the axis is regular (C* = 1), the external conicity parameter C ^ 1 
in general. 



6. Cylinders of incompressible fluid: analytic approach and numerical 
results 

The set of the ordinary differential (field) equations (2.6)-(2.9) can be simplified by 
converting it into the dimensionless form: 

= <"> 

* e=(x + 3 n>§ + ^, (6 ,) 

d 2 „ KN 2 UK 

ds 2 N 2 N 2 y ' 

Here the dimensionless radial coordinate s = (3r, where (3 2 = AttGX^, with fi being 
a characteristic (e.g., the central) fluid density; the dimensionless fluid density is 
X = m/Mo an d the dimensionless pressure II = Xp/fiQ (this is not an EOS). The 
dimensionless metric functions are given in terms of the original metric functions 
U, B,C,F = BC as N = e xu , K = f3B, H = f3F, Q = NH and the dot ' = d/ds. 
When considering an incompressible fluid, we have % = 1, the solution of the equations 
then depends only on the central pressure II C = Xp c /^i , which enters into the solution 
through the initial conditions N(0) = 1, Q(0) = 0, 11(0) = II C , H(0) = 0, H(0) = 1, 
K{0) = 0, K(0) = 1. From now on, we assume C* = 1 on the axis, i.e., we consider 
cylinders without a cosmic string along the symmetry axis. These initial conditions 
lead to a solution regular at s — that yields Q/H — > and thus no problems with 
'0/0' arise during the numerical integration of the differential equations (cf equation 
(6.1)). For a rigorous proof of the existence and uniqueness of a smooth solution in a 
neighborhood of the axis, see Appendix A. 

Before turning to the numerical results on the fully rclativistic cylinders of 
incompressible fluid, it is interesting to derive relativistic corrections to the Newtonian 
cylinders of incompressible fluid analytically. For small radii the analytic expressions 
agree well with fully relativistic numerical calculations. 

In the analytic approach we make Taylor expansions in the dimensionless radial 
distance s around the origin (the axis) s — of all functions entering the field equations 
(6.1)-(6.5), in which, for incompressible fluids, we put \ = 1. Combining the results 
we can express the dimensionless pressure as a series in s in which the coefficients are 
uniquely determined by the value of the pressure at the center: 

n( s ) = n c - i (i + 3n c ) (i + n c ) s 2 + 

+ (6i n c + 21) (i + 3 n c ) (i + n c ) s 4 + o( s 6 ) . (6.6) 

Note that this expansion in s can also be understood as an expansion in the parameter 
A = 1/c 2 since s 2 = (4irG ^o)\r 2 . Hence, it also yields relativistic corrections to the 
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corresponding Newtonian expression (given by the first two terms in the following 
series — see Appendix B): 

p(r) =p c - nGulr 2 + 7rG> r 2 (-7rG^r 2 - 4p)A + 

+ nGr\-^n 2 G^y + l -§^y P 3p 2 )A 2 + 0(A 3 ). (6.7) 
The expansion (6.6) approximates well the numerical results for s <C 1. Assuming 




Figure 1. Pressure fall-off as a function of the distance from the axis. The 
individual curves from upper right to lower left correspond to the values of the 
central pressure El c = 0.001,0.01,0.1,1,10 and 100. The vertical axis gives the 
ratio of the pressure at a given point to the axis pressure U/Tl c . The horizontal 
axis gives the ratio of the coordinate distance r from the axis to the Newtonian 
radius of a cylinder with the same axis pressure: r/J? Nowtonian = r/ p c /wGfj.Q. 
The central pressure increases from upper right to lower left. 

that the dimcnsionless central pressure n c < 1, we can invert the expansion (6.6) to 
find the location of the cylinder's surface, i.e., such value s — 5 at which 11(5) = 0. 
Notice that n c = \p c /fi <C 1 is indeed valid for not very relativistic pressures. 

There are two physical (geometrical) radii of a relativistic cylinder — the proper 
radius R p and the circumferential radius R c , defined by 

f R 1 f s ds 

R ' = I ^ dT = 7Smih N' <6 ' 8) 

" . (6.9) 



Using the expansions of the metric functions K, N in s 2 and substituting for S from the 
condition 11(5) = determined by inverting equation (6.6), we arrive at the following 
results: 

^^l 1 -!^^^). ««»» 
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R c = Jz?h ( 1 -- n c + TT^ n c + o(n3)). (fU1) 



■vTf 17 6415 
7rG/i§ V 8 c+ 1152 

Here the dimensionless central pressure plays the role of an expansion parameter; 
since it involves A, the second and subsequent terms in (6.10), (6.11) determine the 
relativistic corrections. The first terms, of course, coincide: R p = R c = ^pc/irGfi 2 , is 
the radius of a Newtonian cylinder of incompressible fluid with density /iq and central 
pressure p c (see Appendix B). 

Next, defining mass per unit coordinate and proper lengths of the cylinder, Mi 
and M p , by 

GXMi = 2ttGX f nBCe- 3XU dr = \ ( -^ds, (6.12) 

and 

GXMp = 2ttG\ { pBe- 2XU dr = \ [ ^-rds, (6.13) 
respectively, we find expansions 



and 



1 ^ QQ^ 

GXMi = n c - ^n* + ^nl + o(n% (e.w) 



1 Q 7Q1 

GXM P = U C --U 2 C + — Ul + 0(Ui). (6.15) 



Again, the first term yields the expected Newtonian result: Mi — M p — W C /GX 
Pc/Gno- The Levi-Civita mass parameter turns out to be 



m=y§— ^ =2n c (i-In c + ^ + o ( n3)), (6 .i6) 

s=S 

so that in the lowest order GXMi = GXM p = m/2, as is also the case for infinitely 
thin cylindrical shells [6] . 

The conicity parameter, C, outside the cylinder is also uniquely determined by 
the central pressure as it follows from equation (5.5) with C* = 1: 

c = i + K ( 2 m I- - 1) -n3 ( 31n |_ _ |) + o ( n4). (fU7) 

The conicity parameter starts to deviate from 1 only by terms of order 0(11;? lnll c ). 

Let us now turn to the numerical results for fully relativistic cylinders of 
incompressible fluid. The numerical integration of the system (6.1)-(6.5) (with \ = 1) 
enables one to find first the radial distribution of the pressure and then the proper 
and circumferential radii of the cylinders, as well as other physical quantities, all being 
determined uniquely by the central dimensionless pressure II C = p c / floe 2 . 

The dependence of the pressure on the dimensionless distance s from the axis is 
illustrated in figure 1. The higher the central pressure, the more the resulting curve 
deviates from the Newtonian case. The Newtonian curve is represented by a parabola 
y = 1 — x 2 for any central pressure. For low central pressures, the curve is very close to 
this parabola; as the central pressure increases, it drops faster than in the Newtonian 
case and the surface is thus closer to the axis. 

Perhaps the most interesting result can be seen in figure 2. Here the dimensionless 
proper and circumferential radii, R p and R Cl of the cylinders are illustrated for 
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Figure 2. The dependence of the dimensionless proper radius of the cylinder 

R p = GXfioRp and the dimensionless circumferential radius R c = \J GX/ioRc 
of an incompressible fluid cylinder on the dimensionless central pressure II C = 
Xpc/fig- The circumferential radius R c has only a bounded range of values. 



cylinders parameterized by the central pressure Il c . Curiously enough, while the 
proper radius increases with increasing Il c , the circumferential radius starts to decrease 
for higher central pressures, though with still 'physical' values of II C < 1; the maximum 
of the curve i? c (II c ) gives R c ~ 0.213243 and occurs at Il c 0.8 (see figure 2). For 
any finite value of (io = constant > and any finite II C > 0, the coordinate radius of 
the cylinder, R — and, correspondingly, also R p and R c — is finite in accordance with 
the theorem in section 3. 

As discussed below (5.5), knowing the cylinder radius, we can use the matching 
conditions to calculate the Levi-Civita mass parameter, to, and the conicity parameter, 
C, of the vacuum spacetime outside the cylinder for any given Il c . The resulting 
curves are illustrated in figure 3. The Levi-Civita mass parameter, characterizing the 
curvature of the vacuum spacetimes (see equation (2.21)), increases from its flat-space 
value, reaching the magnitude of to m .69 for p c = ^c 2 , and approaching to = 1 for 
the extreme central pressures 1I C 3> 1. Interestingly, the cylinders with still relatively 
low central pressures are so relativistic that they produce Levi-Civita solutions with 
m > 1/2, in which there are no circular timelike geodesies (cf also [20] below Equation 
14) . For an analogous phenomenon in case of static cylindrical shells and their Levi- 
Civita fields, see [6] . 

In the physical region of the pressures, p c < fioc 2 , one can find, by numerical 
interpolation, the following nice analytical approximations (all with a relative accuracy 
better than 10~ 3 up to II C — 1) to the numerical curves illustrated in figures 2, 3: 




(6.18) 
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473 + 2111, 



m 



Pc 



2n c (1 - n c 



i-n c 



225 + 567IT 
733 - 24n c 
413 + 670IL 



c ~ i + n; 



, 4i - 3n c 

: 5 + 22011, 







T 2 168+15n e +5ng 
: 42 + 68n c + lS4n2 



15 
(6.19) 
(6.20) 
(6.21) 




Figure 3. The external Levi-Civita parameter, m, and the external conicity 
parameter, ln(C), as functions of the dimcnsionlcss central pressure, il c = \p c /fio, 
for incompressible fluid (m can attain only values within [0,1)). We set (3 = 
47rGA/io = 1. 



It is not obvious what expression to use as the unit-length mass of the cylinders. 
There are several choices. One can use the Vishveshwara-Winicour definition 
[21] employing Killing vector fields outside the cylinders (in the outer Levi-Civita 
spacetime) and yielding 

m 

M vw = -. (6.22) 

We find Myw £ [0, 1/2). This value crosses the 1/4 limit [6] already for the central 
pressure IT C well below 1. We can also use the Tolman mass 



m T - 2tt^ (/x + X3p)BCe- 2XU dr = ^ U'BC\ r=R - 



2GA N 



, (6.23) 



which is not bounded from above — solutions with unbounded ttit have also been 
found analytically [18]. We can use mass per unit coordinate length, Mi (6.12), with 
no upper bound (see figure 4). There are two more expressions that do exhibit a 
limited interval of values: mass per unit proper length of the cylinder, M p (6.13), with 
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Figure 4. Dimcnsionlcss masses per unit coordinate and unit proper lengths of 
the cylinders, Mi = G\M\ and M v = G\M P , as functions of the dimensionless 
central pressure Tl c . The dimensionless mass per unit proper length M p has an 
upper bound of 1/4. 



GXMp < 1/4 as shown in [22] (see figure 4), and Thome's C-energy scalar defined by 
using the symmetries of the spacetime [23] as 



U= - 



A^A* 

4- 2 iii 2 



1 



fKN 2 d 
\ H ~ds 



H 

N 2 



(6.24) 



where A = 2-xBCe 2U is the area of the cylindrical belt given by r = constant, 
z e [0, 1], <f e [0, 2tt). We find U < 1/8— in accordance with [23] (see figure 5). 

There is a fundamental difference between the two bounded expressions — there 
is no increase in M p outside the solid cylinders but there is non-zero C-energy scalar 
associated also with the outer vacuum Levi-Civita spacetime. The total C-energy 
contained within a cylinder of radius greater than the radius of the solid cylinder is 
given simply by the corresponding expression for the pure Levi-Civita spacetime (no 
integration) U = (1/8) (1 - (1 - m) 4 /r 2 ™ 2 C 2 ) < 1/8. 

We can construct the following invariant expression characterizing the conicity of 
the metric 

,V 'g vv {r\)g r r( 



^ lim Mr*) - R P (n) 
r 2 ^rx R c (r 2 ) - R c {ri) 



= 2- 



-(ri) 



dr 



9w( r i) 



(6.25) 



with R c being the circumferential and R p the proper radius. In Minkowski spacetime 
we find tp = 1. For a cosmic string ds 2 = —dt 2 + dr 2 + dz 2 + ■^■r 2 dip 2 , and we 
have ip = C; in a Levi-Civita spacetime we calculate ip = Cr™ 2 /(l — m) (and thus 
U = (1/8)(1 - [(1 - m)/iP} 2 )). For full cylinders, we obtain 

1 N 



B> - BU> NK - KN 



(6.26) 



-see figure 6. 




Figure 6. The conicity characterizing quantity if) (see equation (6.25)) evaluated 
at the surface of the cylinder with the central pressure n c . 



If we embed a 2-dimcnsional hypcrplane t, z = constant with metric 

ds 2 = c- 2XU (dr 2 + B 2 d^ 2 ) (6.27) 
into a flat Euclidian space with metric 

ds 2 = dTZ 2 + TZ 2 d(p 2 + d( 2 , (6.28) 
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we get 

H = Be-™ = j-. (6.29) 
Equation (6.9) implies 1Z = R c and finally (using (6.8)) 



d^ = 1 i-J dR P i (6 30) 

dR c y (n^k/N)) 2 V dR c 

We conclude 

dR 2 p = dR 2 c + dQ 2 . (6.31) 

Thus i? p measures the length of the embedding curve, see figure 7. It can be shown 
that the embedding surface never degenerates into a cylinder (d(/dR c = oo) at a finite 
distance from the axis. For a cosmic string, this simplifies to a cone d(/dR c = ^/C 2 — 1, 




Figure 7. Embedding of the z,t = constant surfaces according to formula (6.30) 
for central pressures n c = 0.2,1,3,10,25,100, 1000, 10 4 , 10 6 and 10 s , bottom 
to top. The dots on the graph indicate the position of the cylinder's surface. 
R c = \J GXfM)R c is the dimensionless circumferential radius, ( = GA/iqC- 

as expected. Further we find 

Another interesting graph is a plot of unit-length mass within the cylinder as 
a function of the cylinder radius. We have four options: mass per unit coordinate 
or proper length and proper and circumferential radii. In figure 8 we present all 
four quantities. The graph of the unit proper length mass as a function of the 
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circumferential radius resembles the plot of equilibrium spherical configurations where 
we plot the Schwarzschild mass M of the system as a function of its coordinate radius 
R. On the other hand, if we integrate the structure equations in case of a spherically 
symmetric static star of constant density fi , we find 



GXM p = - 



Rn 



— sin AR P 



cos AR, 



^ Q arcsin(.4i? c ) - R c yj\ - A 2 R^j 



(6.33) 



where A — y/WrrGX^Ji, M p is the total proper mass of the star and R p , R c are its 
proper and circumferential radii, respectively (see figure 9). This function is similar to 
the first graph in figure 8 — to the mass per unit coordinate length M\ as a function of 
the dimensionless proper radius R p of the cylinder. There is, however, a fundamental 
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Figure 8. Dimensionless masscs_of the cylinders per unit coordinate and unit 
proper lengths, M\ = G\M\ and M p = GXM P , as functions of their dimensionless 

proper and circumferential radii, R p = yj GXfioRp and R c = yj GXfioR c . The 
graph uses logarithmic scales to reveal the asymptotic behavior of the masses. 

difference between the spherical and cylindrical configurations — in the spherical case 
in standard Schwarzschild spherical coordinates, the radial component of the metric 
tensor depends on the density but not on the pressure. Therefore, we do not need to 
find the pressure as a function of the distance from the center to evaluate M p . For 
spheres of incompressible fluid, the dependence of the proper mass on the distance 
from the center inside the sphere is common for all spheres and it is the same as 
the dependence of the total proper mass on the radius of the spheres. This is not 
the case for cylinders of incompressible fluid where each value of the central pressure 
determines a unique curve p(r). 
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Figure 9. Dimensionless total proper mass M v = yj GA/io GXM P of a static 
spherical star of constant density /xq as a function of its dimensionless proper 
radius R p = R p and dimensionless circumferential radius R c = 

y^7rGA/irj/3 R c ■ In the spherical case, the circumferential radius R c coincides 
with the coordinate radius R in the standard spherical coordinates. 



7. Concluding remarks 

Although the main results were announced in the introductory section already, let us 
briefly summarize them here, emphasizing those aspects that appear to be new in the 
subject of cylindrical symmetry in general relativity. Our primary goal, in contrast to 
majority of other works, has been (i) to understand the global character of spacetimes 
with static perfect fluid cylinders, (ii) to study weakly gravitating cylinders and their 
Newtonian limit and (iii) to analyze in detail cylinders of incompressible fluid. Except 
for the last item, we did not start from an a priori equation of state, so our results 
are of a general character. 

We have shown that for any smooth monotonic equation of state and for any 
density at the axis of symmetry there exists a unique solution in some neighborhood 
of the axis which is regular at the axis and can be uniquely extended to a global 
solution. In other words, the equation of state determines a one-parameter family of 
global spacetimes. If the pressure vanishes at a finite value of the radial coordinate 
then a unique Levi-Civita vacuum solution can be joined to the inner perfect-fluid 
solution. In particular, we prove that the cylinder has a finite radius if the monotonic 
equation of state admits a non-vanishing density for zero pressure (as is the case 
for incompressible fluid, for example). In general, the outside Levi-Civita solution is 
determined by both the mass (curvature) parameter and the conicity parameter. The 
need for a nontrivial conicity parameter in the external vacuum region and the fact 
that both the mass and conicity parameters are determined uniquely by the value 
of the density /pressure at the axis of the fluid cylinder do not appear to have been 
elucidated before. 

In the Newtonian limit, we prove that the spatial metric inside the fluid cylinders 
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is conformally flat. The relativistic Tolman mass becomes the Newtonian mass per 
unit length. For comparison, we also discussed cylinders and spheres of perfect fluid 
constructed ah initio in the Newtonian theory. 

In the case of relativistic cylinders of incompressible fluid, no analytical solution 
is available. However, we succeeded in deriving analytic results for the relativistic 
corrections to the Newtonian cylinders. For weekly relativistic cylinders, the Levi- 
Civita mass parameter and the conicity parameter, for example, are given by equations 
(6.16) and (6.17). For fully relativistic cylinders, various physical quantities of 
interest were found by numerical integration of the field equations and they are 
exhibited graphically (figures 1-8 in section 6). The resulting configurations are 
determined uniquely by the dimcnsionless pressure Il c = p c / floe 2 at the axis. As 
the central pressure increases, the pressure away from the axis decreases more rapidly 
so that relativistic cylinders become more compact than the Newtonian ones. A 
remarkable phenomenon, unnoticed so far, arises in the relativistic regime: while 
the (dimcnsionless) proper radius of the cylinders increases with increasing central 
pressure, the (dimensionless) circumferential radius starts to decrease for still 'physical' 
values of II C < 1; the circumferential radius has only a bounded range of values. 
This phenomenon is nicely illustrated by the embedding diagrams of the surfaces 
z, t ^constant (figure 7). The external Levi-Civita mass parameter m increases with 
II C , approaching m — 1 as II C 3> 1, whereas the quantity tp characterising the conicity 
increases without limit as IT C — > oo. The dimensionless mass per unit proper length 
inside the cylinders starts to decrease for highly relativistic cylinders (it has an upper 
bound of 1/4), but the analogous mass per unit coordinate length increases without 
bound. Thome's C-energy increases with n c and approaches 1/8 as n c becomes large. 

At the end of the preceding section, a comparison of cylindrical and spherical 
configurations of incompressible fluid shows similarities and fundamental differences 
between the two cases (figures 8 and 9). 

The 'external' Levi-Civita conicity parameter, the importance of which has been 
emphasized throughout the text, takes various values depending on the 'inner' cylinder, 
which is regular everywhere. However, to include more general situations, we also 
considered an infinitely thin cosmic string (see equation (2.2)), i.e., a conical singularity 
along the axis of symmetry inside the cylinder. This produces the 'axis conicity' which 
enters various formulae, like the Tolman mass, but it does not, for example, influence 
the conclusions of the theorem stating that a cylinder has a finite radius if the equation 
of state admits positive density at vanishing pressure. 

Finally, to 'end with the beginning', let us remark again that only full 
understanding of the static situation will enable a sufficiently thorough treatment 
of problems such as the interaction of cylindrical gravitational waves with static 
cylindrical matter configurations. 
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Appendix A. Regularity of the axis 

In order to prove the existence of solutions regular on the axis we begin with the 
following two equations (2 x (2.7) + B x (2.8) + C x (2.9)) and (2.6) written for 
F = BC: 

F" = &Ti\ 2 Ge- 2UX ApF = V(U, X)F , (A.l) 

U" + ^U' = AirG(p + X 2 3p)e- 2UX = W(U, A) . (A.2) 
r 

For a given EOS, /x(p), and the value of the pressure on the axis, p c > 0, the function 
p(U) is uniquely determined from equation (2.10) since p + n does not change the sign. 
The boundary values at r = arc F(0) = U(0) = U'(0) = 0, F'(0) = 1. 
We want to write the above system in the form 

x^ + Yf = xG(x,f(x))+g(x), (A.3) 

where Y is a constant n x n matrix and f(x) a vector. For this form, the existence and 
uniqueness of a smooth solution in a neighborhood of x = was shown in theorem 1 
in [24], provided the matrix Y has positive eigenvalues. 
First we define 

F(r)=rf(r 2 ), (A.4) 
where /(0) = 1 and /'(0) is finite, and obtain the equation 

U" + -U' + 2rC = W. (A.5) 

r f 

With x — r 2 as radial coordinate and U(r 2 ) = U(r), we get 

4xU" + 4U' + 4xLtj' = W(U,X), (A.6) 

where ' is the derivative with respect to the argument of the function, i.e, x. We write 
this equation as a first order system 

U' = v, (A.7) 
f' 

Axv' + Av + 4xjv = W(U, A) . (A.8) 

With U(x) = xu(x), we obtain 

xu' + u-v = Q. (A. 9) 
Similarly, we rewrite the equations for / as a singular first-order system. Defining 

f' = 9, (A. 10) 

we get 

Axg + &g = Vf . (A. 11) 
To make the equation singular, we define 

f = xh+l (A. 12) 
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and obtain 

xh' + h-g = 0. (A.13) 
Now we finally have the system in the desired form of equation (A. 3): 

xbl + h 1 

xi)' + v = -x— v + -W(xu, A) , (A. 14) 

xh + 1 4 

xu' + u-v = 0, (A.15) 
xg' + \g = j(l + xh)V(xu,\), (A. 16) 

(A.17) 



(A.18) 



which has positive eigenvalues 1, |. 

Appendix B. Static Newtonian configurations of matter 

For a Newtonian cylinder, the balance equation reads 

p' + g» = 0, (B.i) 

where g is gravitational acceleration, p is pressure and /i is the fluid density (variable 
in general). Using the Gauss theorem, we have 

pl + ^GH j\rdr = Q. (B.2) 
r Jo 

Since there appears an integral quantity, namely the total mass contained within the 
cylinder up to radius r, we obtain a second-order differential equation 

AnGrfi 3 = rp'p! - fj,p' - p,rp" . (B.3) 

For incompressible fluid of density /x , pressure p is a quadratic function of distance 
from the axis and, therefore, the mass Mi per unit length of the cylinder and its radius 
R are given by 

p = p c - nG^r 2 , 



Pc 

TrG/xg' 



R = (B-4) 



M 1= Pc 



Is it possible to find the dependancc of the cylinder radius on the central pressure 
in other cases as well? Let us consider, for example, a polytropic EOS, p — ct[i K , 
with a dimensionless polytropic index n. Substituting the EOS into equation (B.3), 
we have 

4irGrfi 3 = anr(2 - k)^ 1 ^' 2 - ok/ p! - anrp^p!'. (B.5) 
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One of the solutions is fi — ( — a ^^_^ — ^ • This, however, is either zero (k > 2) on 
the axis or divergent (k < 2) and thus it is not physically relevant. For k = 2 we find 



the non-divergent solution to be the Bessel function of the first kind Jo(y ^r" r )- I n 
this case the radius of the cylinder does not depend on the central pressure. 

This peculiar result can be deduced using dimensional analysis as follows. The 
radius R of a Newtonian polytropic self-gravitating cylinder is described completely 
by three quantities: p c , /i c and G. From these quantities only one dimcnsionless 
parameter — it = GR 2 fi 2 p~ x — can be constructed. Thus, the Buckingham ir theorem 
[25] says, that the physical law must take the form /(k,7t) = 0, and, assuming 
uniqueness, we can write 7r = g 2 (n), i.e. 



Pc n -\ ( \ P 

Since Pc/f-c is constant for a given EOS, the derivative 
dR k-2 R 



R = ffW t /^ = G-*g{K)[ ^)p—. (B.6) 



(B.7) 



dp c 2k p c 

suggests that if a unique solution of hydrostatic equations exists, then the sign of 
dR/dp c will be the same as the sign of k — 2. Moreover, for k = 2 the radius of the 
cylinder will not depend on p c at all. 

We note that the same argument cannot be used in the relativistic case since 
then we have another quantity, namely c, the velocity of light, which we need to 
take into account. In addition to the dimcnsionless quantity ir mentioned above, we 
can also construct the dimensionlcss central pressure n c = p c jc 2 [i c . Consequently, 
/(k,7t, n c ) = 0, and we do not find the derivative (B.7) explicitly. 

To demonstrate the dependance of these results on the symmetry of the 
configuration, we give the results for Newtonian spheres. The balance equation reads 

AirGrfi 3 = rp'fi' - 2fip - ^rp" . (B.8) 

For incompressible fluid of density /j,o, pressure is again a quadratic function of the 
distance from the center. The total mass M contained within the sphere and its radius 
R read 

P =Pc-^Gn 2 r 2 , 

I- 77^3 (B.<>) 




The dimensional analysis for a polytropic EOS yields exactly the same results as in 
the cylindrical case. 



Appendix C. Conicity outside general perfect-fluid cylinders 

We have seen that the conditions at the axis determine uniquely the fields both inside 
and outside the cylinders. In particular, we can choose the axis to be regular or to 
have an inner conicity C* ^ 1, and the conicity parameter C of the external field is 
determined uniquely (cf equation (5.5)). In section 6, we have shown numerically that 
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the conicity parameter C outside relativistic self-gravitating cylinders of incompressible 
fluid is always greater than 1. Is that the case for any cylinder of perfect fluid? Let 
us investigate what happens if one just puts C = 1 in the exterior Levi-Civita metric. 
We consider the cylinders discussed by Stela and Kramer [9] and use their coordinate 
system. They start out from the following metric describing a solid cylinder composed 
of perfect fluid: 

ds 2 = _ e 2x dt 2_ + y_l^l dx 2 + e -2x ,2/^2 + ^ 2 N 

8irp v ' 

where t_, x, £ and y>_ G [0, 2n) are the temporal, radial, axial and azimuthal 
coordinates, respectively, p is pressure, and y, z, k, h are functions of x. This applies 
to the range x € [0, x\), where x\ denotes the surface of the cylinder, p{x\) = 0. 
Outside, they use Levi-Civita spacetime in the form 

ds 2 = -p 2m dt\ + p- 2m \p 2m2 (dz 2 + dp 2 ) + p 2 dip 2 + ] , (C.2) 

i.e., without a general conicity parameter, so the metric (C.2) is our metric (2.24) with 
C = 1. It is claimed in [9] that these two spacetimes, described by (C.l) and (C.2), 
can be smoothly joined on the surface p = 0. However, this cannot be achieved if the 
complete spacetime is to represent a solid cylinder and its external gravitational field 
with a fixed range of the angular coordinate, ip, throughout the whole spacetime. We 
have to consider a general conicity parameter C ^ 1. 

To demonstrate this, we shall use Israel's formalism [19]. We consider two 
hypersurfaces located at x = x\ inside and at p = p + in the outside spacetime. 
These two hypersurfaces are intrinsically flat and, therefore, we may identify them 
locally. Let the coordinates on the resulting hypersurface (representing possibly a 
matter shell) be chosen as T G (—00,00), Z e (—00,00) and $ <G [0, 27rp^_ _m ) so that 
the induced 3- metric is flat: 

t_e Xl =T = t+p™, 

£ e k( Xl )- Xl =Z = zpf m - 1 \ (C.3) 

The induced surface (3-dimensional) energy-momentum tensor calculated from the 
jump of the extrinsic curvatures (see [19]) then turns out to be 



SttStt = ;/-^r (y + z-2)- p™- m2 -\l - mf 
yz 1 



8nb$$=p + m - W y, 

V yz - 1 

where functions y and z are evaluated on the surface of the cylinder (p = 0). Einstein's 
equations yield the following relations for the metric functions 

y = (l-yz){Fy-2), z = (1 - yz){Tz - 2), 

= P +(p + p), 

where T — (p + 3p)/2p (see equations (1), (2) in [9]). On the surface, G xx — then 
implies yz = 1, and the fraction under the square root symbols in (C.4) is thus of the 
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'0/0' type. Using (C.5) and l'Hospital's rule, we find on the surface 

z = y/y 2 , (C.6) 

8irp / Airpy 



yz — 1 V y 

If the inner solution is to be smoothly extended into the Levi-Civita spacetime without 
any matter shell on the surface (no 'surface layer'), Sab must be zero. Combining 
equations (C.4) and (C.6), we get on the surface 



4 -^(y+ 1 --2) = P r m2 - 1 (l-mf, 

y y 

= P ™- m '-\ (C.7) 



Anny 1 ~_ m *_i 



y y 



m-m 2 -l 2 

y = p+ rn . 

y 

Comparing the last two equations, we see that y{x\) = ±m = l/z{x\). Then the first 
equation yields 

y(x 1 )=m = l/z(x 1 ). (C.8) 
This leaves us with only one equation fixing the outer radius of the cylinder: 



If one is only interested in the value of the Levi-Civita parameter m this equation is 
sufficient. However, there is one more (global) junction condition, namely that the 
proper length of a hoop with constant T, Z must be the same as measured on both 
sides of the surface of the cylinder: 

2ne h( Xl )- Xl = 2 7rp 1 + - m , (CIO) 

so 

h{x l ) = x 1 + (1 - m)lnp + . (Cll) 

Integrating Einstein's equations from the axis to the surface where p{x\) — 0, we 
infer the value of m outside just by reading off the value of y{x\) — 1/ z(x\) as seen 
from equation (C.8). However, to ensure that there be no surface layer at p + , wc 
have to require that also conditions (C.9) and (Cll) be fulfilled. There are thus two 
conditions but only one free parameter — the outer radius p + . Combining equations 
(C.9) and (Cll), we see that any solution must satisfy the relation 



.^ = V ye^)^, (C.12) 

y 

where all quantities are evaluated at the surface x = x\. For a given EOS, the values of 
all these quantities are determined by a single parameter — the central pressure. This 
means that, generally, we cannot choose the central pressure arbitrarily since it has 
to satisfy (C.12). Therefore, only special cylinders will comply with this condition. 
Moreover, we have shown in section 6 that in case of incompressible fluid one gets 
C > 1 for any central pressure (see figure 3 and expansion (6.17)) and thus none of 
the cylinders composed of incompressible fluid can satisfy (C.12) because we derived 
this equation starting from the Levi-Civita metric (C.2) with C = 1. 
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Surprisingly enough, there exist EOSs for which special cylinders can be 
constructed that do satisfy this condition and hence the outside Levi-Civita metric 
with C = 1 applies. To give an explicit example, let us consider the Evans analytical 
solution [8] , which is also discussed in [9] . The equation of state is n = /io + 5p; /io > 
is (arbitrarily chosen) density at the surface where p = 0. In the coordinate system of 
[9], the solution can be written as follows: 

P( X ) - TT, ..,_M" „-(te _ X 



6 V 4 

,2 









^ 4 


-0- 




- 4e 3x 





Mo 

n c = \ ( ^ - 1 ) , e 3 - = vTTmt 



1 - 4e 3:c 

^ 2( fl2 - 4e3^) ' = 2(1-4^) ' (C ' 13) 

, , 1 2y/l + 6II C - 2 1 
2/(^1 ) = -7 — 7 = m = , 7l , r < 



lim 



z{x 1 ) 4-y/l + 6n c - 1 2' 

- 1 1 18a 



+zi p /i (2a — l)(a — 2) ' 

a is an arbitrary constant determining the central dimensionless pressure. The 
condition (C.12) gives a relation 6 between the parameter a (or, equivalently, II C ) and 
the surface density /io that has to be satisfied in order that C = 1: 

[ M o/(a - 2)} {a - 2)2 = 2 2 ^ 2+2a - 2 h- 7a2+a - l a- 3a2 [TT{2a - i)3](2a-i)(a+i). (ai4) 

Therefore, given a specific equation of state, which fixes the value of /x , the pressure on 
the axis cannot be chosen arbitrarily if the outside field has C = 1. Not all physically 
admissible cylinders can be sources of the Levi-Civita field with C = 1. Only if p c (i.e. 
n c = Pc/mo) is chosen so that relation (C.14) is indeed satisfied, the external field has 
C = 1. However, any value a > 2 gives a positive value of both the central pressure, 
II C , and surface density, /io, and it thus gives a cylinder with C = 1 outside. 
If we admit C ^ 1 outside, (C.ll) changes to 

h{x\) = X\ + (1 — to) In p + — InC, (C.15) 

and we can smoothly join any cylinder composed of perfect fluid to the outer 
Levi-Civita solution. On the other hand, most characteristic quantities are either 
independent of /xo (e.g., the Levi-Civita parameter to, mass per unit coordinate or 
proper length of the cylinder), or they scale with /zo (R p = r p (IL c ) / \fjMy) so that it 
is possible to use a single cylinder with a given II C for the description of properties 
of a whole family of cylinders. Nevertheless, one does not determine the complete 
spacetime geometry outside the cylinders without knowing C. If, for example, one 
studies the focusing of light rays passing on both sides of a fluid cylinder, one needs 
to know C. 



References 

[1] Bicak J 2000 Selected solutions of Einstein's field equations: their role in general relativity and 
astrophysics Einstein's Field Equations and Their Physical Implications. Selected Essays 
in Honour of Jiirgen Ehlers (Lecture Notes in Physics vol 540 ) ed B G Schmidt (Berlin: 
Springer) pp f-f26 

6 The complicated nature of expression (C.14) follows from the fact that the requirement of C = 1 
outside the cylinder is rather artificial. 



Static fluid cylinders and their fields: global solutions 



28 



Bondi H 2000 Proc. R. Soc. A 456 2645 

Mena Marugan G A 2000 Phys. Rev. D 63 024005 

Manojlovic N and Mena Marugan G A 2001 Class. Quantum Grav. 18 2065 
Goncalves S M C V 2002 Phys. Rev. D 65 084045 

Goncalvcs S M C V 2003 Class. Quantum Grav. 20 37 
Bonnor W B 1999 The static cylinder in general relativity On Einstein's Path: Essays in Honor 

of Engelbert Schucking ed A Harvey (New York: Springer) pp 113-9 
Bicak J and Zofka M 2002 Class. Quantum Grav. 19 3653 

Stephani H, Kramer D, Maccallum M A H, Hocnsclaers C and Herlt E 2003 Exact Solutions to 

Einstein's Field Equations 2nd edn (Cambridge: Cambridge University Press) 
Evans A B 1977 J. Phys. A: Math. Gen. 10 1303 
Stela J and Kramer D 1990 Acta Phys. Pol. B 21 843 

Carot J, Senovilla J M M and Vera R 1999 Class. Quantum Grav. 16 3025 
Kramer D 1987 Class. Quantum Grav. 5 393 

Teixeira A F, Wolk I and Som M M 1977 Nuovo Cimento B 41 387 

Bronnikov K A 1979 J. Phys. A: Math. Gen. 12 201 

Haggag S and Desokey F 1996 Class. Quantum Grav. 16 3221 

Scheel M A, Shapiro S L and Tcukolsky S A 1993 Phys. Rev. D 48 592 

Ehlers J 1997 Class. Quantum Grav. 14 A119 

Ehlcrs J 1998 The Newtonian limit of general relativity Understanding Physics ed A K Richtcr 

(Katlenburg, Lindau: Copernicus Gesellschaft e.V.) pp 1-13 
Philbin T G 1996 Class. Quantum Grav. 13 1217 
Israel W 1966 Nuovo Cimento B 44 1 

Israel W 1967 Nuovo Cimento B 48 463 (erratum) 
Bonnor W B 1991 Gen. Rel. Grav. 24 551 

Vishveswara C V and Winicour J 1977 J. Math. Phys. 18 1280 
Anderson M R 1999 Class. Quantum Grav. 16 2845 

Thornc K S 1965 PhD Thesis Princeton University, available from University Microfilms Inc., 
Ann Arbor, MI 

Rcndall A D and Schmidt B G 1991 Class. Quantum Grav. 8 985 
Buckingham E 1914 Phys. Rev. 4 345 



